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Abstract 

Inverse modeling of nitrogen oxide (NO x ) emissions using satellite-based N0 2 obser- 
vations has become more prevalent in recent years, but has rarely been applied to 
regulatory modeling at regional scales. In this study, OMI satellite observations of N0 2 
5 column densities are used to conduct inverse modeling of NO x emission inventories for 
two Texas State Implementation Plan (SIP) modeling episodes. Addition of lightning, 
aircraft, and soil NO x emissions to the regulatory inventory narrowed but did not close 
the gap between modeled and satellite observed N0 2 over rural regions. Satellite- 
based top-down emission inventories are created with the regional Comprehensive Air 
10 Quality Model with extensions (CAMx) using two techniques: the direct scaling method 
and discrete Kalman filter (DKF) with Decoupled Direct Method (DDM) sensitivity anal- 
ysis. The simulations with satellite-inverted inventories are compared to the modeling 
results using the a priori inventory as well as an inventory created by a ground-level 
N0 2 based DKF inversion. The DKF inversions yield conflicting results: the satellite- 
15 based inversion scales up the a priori NO x emissions in most regions by factors of 1 .02 
to 1.84, leading to 3-55% increase in modeled N0 2 column densities and 1-7ppb 
increase in ground 8 h ozone concentrations, while the ground-based inversion indi- 
cates the a priori NO x emissions should be scaled by factors of 0.34 to 0.57 in each 
region. However, none of the inversions improve the model performance in simulating 
20 aircraft-observed N0 2 or ground-level ozone (0 3 ) concentrations. 


1 Introduction 

Nitrogen oxides (NO x = NO + N0 2 ) in the troposphere are primary air pollutants, emit- 
ted from both anthropogenic sources like fossil-fuel combustion and biomass burning, 
and natural sources such as soil microbial processes and lightning. NO x also acts as 
25 a precursor of a secondary air pollutant, tropospheric 0 3 , when it reacts with the oxi- 
dation products of volatile organic compounds (VOCs) in the presence of sunlight. Ox- 
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idation with hydroxyl (OH) radical is the dominant sink of NO x , leading to atmospheric 
nitric acid (HN0 3 ) formation. The atmospheric lifetime of tropospheric NO x varies from 
a few hours in summer to a couple of days in winter (Seinfeld and Pandis, 2006). 

NO x emission inventories used in air quality modeling are typically developed by 
5 a bottom-up approach based on estimated activity rates and emission factors for each 
category. Due to inaccuracies in determining these rates and factors, the uncertainty 
in NO x emission inventories has been suggested to be as high as a factor of two and 
classified as one of the top uncertainties in ozone simulations and sensitivity analysis 
(Hanna et al., 2001; Xiao et al., 2010). 

io Inverse modeling techniques can be used with atmospheric models to estimate 
model variables that may not be directly measurable (Gilliland and Abbitt, 2001). In- 
verse modeling generates an optimized “top-down” NO x emission inventory for air qual- 
ity models by minimizing the difference between observed and modeled N0 2 concen- 
trations, providing an opportunity to identify possible biases in the bottom-up NO x emis- 
i 5 sion inventory (Napelenok et al., 2008). However, as uncertainties may also associate 
with the measurement data and the inverse methods themselves, inverse modeling has 
its own limitations. Hence, it is valuable to compare both bottom-up and top-down NO x 
emission inventories in order to improve the understanding of NO x emissions. 

Several inverse modeling studies have used surface N0 2 measurements (Mendoza- 
20 Dominguez and Russell, 2000; Quelo et al., 2005; Pison et al., 2007) or aircraft N0 2 
measurements (Brioude et al., 201 1 ) to constrain NO x emissions. Compared to ground 
and aircraft measurements, satellite-based observations generate greater spatial cov- 
erage of N0 2 . Studies on combining satellite N0 2 measurements with inverse modeling 
techniques to create the top-down NO x emission inventories also have been conducted 
25 recently in both global scale (Martin et al., 2003; Muller and Stavrakou, 2005; Jaegle 
et al., 2005; Lin et al., 2010) and regional scale modeling (Konovalov et al., 2006, 
2008; Deguillaume et al., 2007; Napelenok et al., 2008; Kurokawa et al., 2009; Zhao 
and Wang, 2009; Chai et al., 2009). 
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Discrete Kalman filter (DKF) (Prinn, 2000) is an inverse modeling method that solves 
the inverse problem iteratively, and can be applied to the cases with linear or weakly 
non-linear relationships between emissions and pollutants. It has been used in several 
studies to constrain emissions of carbon monoxide (Mu I hoi land and Seinfeld, 1995), 
5 chloroflourocarbons (Haas-Laursen et al., 1996), isoprene (Chang et al., 1996) and 
ammonia (Gilliland et al., 2003). Most recently, Napelenok et al. (2008) applied the 
DKF method to the regional Community Multiscale Air Quality (CMAQ) model, gen- 
erating a top-down NO x emission inventory for the southeastern United States using 
Scanning Imaging Absorption Spectrometer for Atmospheric Chartography (SCIAM- 
io CHY) (Bovensmann et al., 1999) satellite N0 2 data. 

Despite the growing number of scientific studies conducting satellite-based inver- 
sions of NO x emissions, the applicability of these methods to state-level regulatory 
attainment modeling has not been widely explored. In this work, the DKF method in- 
troduced by Napelenok et al. (2008) is applied with finer resolution satellite N0 2 data 
is now available from the Ozone Monitoring Instrument (OMI) as well as ground-level 
N0 2 observations, to constrain NO x emissions for actual regulatory modeling episode 
in Texas. Lightning and aircraft NO x emissions are added to the base case NO x emis- 
sion inventory to address the bias noted by Napelenok et al. (2008) of regional models 
underestimating upper tropospheric NO x . The DKF inverted a posteriori emissions are 
20 compared to the base case emissions, the a priori emissions and a posteriori emissions 
derived by the inversion method of Martin et al. (2003). 


2 Methodology 

2.1 Model inputs and configurations 

Base case model inputs were taken from episodes developed by the Texas Commis- 
25 sion on Environmental Quality (TCEQ) for Texas ozone attainment planning. CAMx ver- 
sion 5.3 (ENVIRON, 2010) was used in this study to simulate two modeling episodes 
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in 2006 with high ozone concentrations in the Dallas-Fort Worth (DFW) region, from 
31 May to 1 July, and in the Houston-Galveston-Brazoria (HGB) region, from 13 Au- 
gust to 1 5 September (Fig. 1 ). The NCAR/Penn State (National Center for Atmospheric 
Research/Pennsylvania State University) Mesoscale Model, Version 5, release 3.7.3 
(MM5v.3.7.3) (Grell et al., 1994), conducted with the Eta-PBL scheme, was used to 
generate the meteorological fields with 43 vertical layers. The preprocessor MMSCAMx 
was used to convert MM5 outputs into CAMx-ready meteorology inputs. The ver- 
tical configuration of CAMx modeling consists of 17 vertical layers for the August- 
September modeling episode, whereas 28 vertical layers were used for the June mod- 
eling episode. Modeling was conducted with the Carbon Bond version 2005 (CB-05) 
chemical mechanism, PPM advection scheme, and K-theory vertical diffusion scheme 
(TCEQ, 201 0, 201 1 ). Boundary conditions for the 36 km eastern US domain were gen- 
erated by the Model for Ozone and Related Chemical Tracers (MOZART) global model 
(ENVIRON, 2008). 

2.2 Emission inventory 

Base case emission inventories were provided by TCEQ. The point source emissions 
were from the State of Texas Air Reporting System (STARS) database which col- 
lects emission information from approximately 2000 point sources annually, and the 
EPA’s acid rain database (ARD) which contains emissions from electric generating 
units (EGUs). The on-road mobile emission inventory was generated by Motor Vehicle 
Emission Simulator 201 0a (MOVES201 0a), and the non-road mobile inventory was de- 
veloped by National Mobile Inventory Model (NMIM) and the Texas NONROAD (TexN) 
mobile source model. The area source inventory was projected by the EPA Economic 
Growth Analysis System (EGAS) model based on 2005 emissions from the Texas Air 
Emissions Repository (TexAER) database. The Emission Processing System, version 
3 (EPS3) (ENVIRON, 2007) was used for processing the point, mobile, and area emis- 
sions to the model-ready format (TCEQ, 2010, 2011). Biogenic emissions were gener- 
ated by the Global Biosphere Emissions and Interactions System (GloBEIS) biogenics 
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emissions model, version 3.1 (Yarwood et al., 1 999), with soil NO x emissions estimated 
by the Yienger and Levy method (Yienger and Levy, 1995). 

Lightning and aircraft NO x emissions in the upper troposphere were missing in the 
base case emission inventories and should be added before conducting inversions. 
In this study, lightning NO emissions were developed based on National Lightning 
Detection Network (NLDN) data obtained from Vaisala Inc., following the approach 
of Kaynak et al. (2008). Intra-cloud lightning flashes were treated as three times of 
cloud-to-ground lightning flashes with 500 mol NO emission per flash. Lightning NO 
was placed into the model to match the time and location of NLDN flashes, and then 
distributed vertically based on the profile obtained from the mean April to Septem- 
ber 2003-2005 vertical distribution of VHF sources from the Northern Alabama Light- 
ning Mapping Array (Koshak et al., 2004). Global aircraft NO x emissions of year 2005 
in 0.1° x 0.1° resolution were obtained from the Emission Database for Global Atmo- 
spheric Research (EDGAR) v4.1 (http://edgar.jrc.ec.europa.eu/datasets_grid_list41. 
php?v=41&edgar_compound=NOx) and mapped to our modeling domain and placed 
at 9 km altitude. 

2.3 Inversion regions 

Five urban areas: Houston-Galveston-Brazoria (HGB), Dallas-Fort Worth (DFW), 
Beaumont-Port Arthur (BPA), Northeast Texas (NE Texas), and Austin and San An- 
tonio; plus two surrounding rural areas: North Rural area (N rural) and South Rural 
area (S rural) (Fig. 1) were designed as inversion regions for the DKF inversions of 
NO x emissions. The five urban regions are all air quality planning areas included in 
Texas SIP development (Gonzales and Williamson, 2011). HGB and DFW were clas- 
sified by US EPA as ozone nonattainment areas for violating the 1997 ozone National 
Ambient Air Quality Standard (NAAQS) of 84ppb. BPA was designated as an ozone 
maintenance area, and NE Texas, Austin and San Antonio were designated as ozone 
early action compact areas under that standard. However, the recent tightening of the 
NAAQS to 75ppb has heightened interest in ozone reduction in all of these regions. 
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The sensitivities of N0 2 concentrations to boundary conditions and to NO x emissions 
from each inversion region and the border region (the area between model boundary 
and inversion regions) were computed through DDM. The border region minimizes the 
impacts from boundary conditions on the inversion regions to the level of only 2%. 
The DDM sensitivities show that NO x emissions from each urban region has the most 
impact on N0 2 concentrations within that region, and has less than 10% influence on 
other regions. 
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2.4 Inversion methods 

Two methods are applied for inverse modeling: a direct scaling method introduced by 
Martin et al. (2003), and the DKF method. 

2.4.1 Direct Scaling (DS) inversion method 

The DS method applies the ratio between satellite N0 2 observations and modeled N0 2 
concentrations to scale the bottom-up NO x emissions in each grid cell: 

Q s 

( 1 ) 

where E x is the top-down NO x emission rate; E b is the bottom-up NO x emission rate; 
Q s and Q m are the satellite and modeled N0 2 column densities, respectively. 

This method was developed in a global model with coarse grid resolution and as- 
sumes that the N0 2 concentration in each model grid will not be affected by the NO x 
emitted from surrounding grids. However, in a regional model with relatively small grid 
size, this assumption may fail, generating a spatial smearing error when NO x lifetime 
is longer than the horizontal transport time (Martin et al., 2003; Boersma et al., 2008; 
Lamsal et al., 2010; Turner et al., 2012). Martin et al. (2003) indicated that the spatial 
smearing error can be neglected if the grid length is greater than 100 km. Therefore, 
smoothing kernels (Toenges-Schuller et al., 2006; Boersma et al., 2008; Lamsal et al., 
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2010) need to be applied in order to alleviate the spatial smearing error in CAMx by 
accounting for the emissions from adjacent grid cells in developing the top-down NO x 
emissions. The smoothing kernel is defined as 


K = 


1 

k + 8 



( 2 ) 
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where k is a smoothing parameter, and is determined by applying the smoothing kernel 
(K) to each grid cell in the bottom-up NO x emission inventory with different k values 
until the correlation between smoothed bottom-up NO x emissions and corresponding 
CAMx modeled N0 2 column density reaches a maximum. The smoothing kernel (K) is 
then applied to Eq. (1) to form Eq. (3), 






£k 



(3) 


where / and j represent column and row in horizontal model grids. 


2.4.2 DKF inversion method 


The direct scaling inversion approach, as described above, creates spatial smearing 
errors when applies to the regional models with fine resolution. Meanwhile, it assumes 
concentrations scale proportionally with emissions; hence, the nonlinearity between 
N0 2 concentrations and NO x emissions becomes problematic because NO x may influ- 
ence its own lifetime by influencing concentrations of OH radicals (Martin et al., 2003). 
The DKF inversion (Fig. 2), however, solves the spatial smearing problem by taking 
the spatial relationship between N0 2 concentrations and NO x emissions directly from 
model simulations, and also reduces the non-linearity issue by performing the inversion 
iteratively. 


CO 

O 

c 

CO 

CO 

o' 

“0 

O 

T3 

CD 


W. Tang et al. 


Title Page 


Abstract ■ Introduction 


Conclusions | References 


co 

o 

c 

CO 

CO 

o' 

o 

"O 

0) 

T3 

CD 


Tables 


Figures 



Close 


Full Screen / Esc 


co 

o 

c 

CO 

CO 

o' 

o 

"0 

a> 

"O 

CD 


Printer-friendly Version 


Interactive Discussion 


(cc) © 


17486 















To constrain NO x emissions, the DKF inversion includes two processes at each 
time step: the measurement update (correction) process and the time update (pre- 
diction) process (Rodgers, 2000; Welch and Bishop, 2001). In the measurement up- 
date process at time step k (Eqs. 4-6), the inversion corrects the predicted NO x emis- 
5 sion (E “ 0 k ) and error covariance (P ~ 0 k ) by incorporating the measurement data 

(C— d ) ancl Kalman Gain (G k ), and then generates the corrected emission (E NOx ^) 

and error covariance (P NOx ,J- 


” ^NO x ,k^k + ^ k ) 

(4) 

c , r* //■'measured /■*modeled\ 

>=no x ,/c - t N0 ^ k + w \y mzik ~ So 2 ,/c ) 

(5) 

10 Pno x ,/c = (■ ~ G/(S^)Pno x ^ 

(6) 


S represents the N0 2 sensitivity to NO x emissions. R is the measurement error co- 
variance, and it relates to the uncertainties in OMI and ground N0 2 measurements. In 
here, the uncertainty for the AQS ground N0 2 measurements was set to 0.15 (US EPA, 
is 2006) and for the NASA standard OMI N0 2 , version 2, was set to 0.3 (Bucsela et al., 
2013) for all diagonal elements in R. The error covariance (P) relates to the uncertainty 
in the NO x emission inventory, and the uncertainty value of 2.0 (Hanna et al., 2001; 
Napelenok et al., 2008) was chosen here for all diagonal elements in P. To simplify, 
off-diagonal elements in R and P were set to zero, because we assume each inversion 
20 region is an independent element. 

In the time update process at time step k, the inversion process predicts the emission 
(E no ^ +1 ) and the error covariance (P ~ 0 k+ ^) for the measurement update process at 

time step k + 1 , based on the corrected emission {E m ^ k ) and error covariance (P NOx ,/() 
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from the measurement update process at time step k (Eqs. 7, 8). 


(7) 

( 8 ) 


ACPD 
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E NO x ,/c+1 “ M /c E NO x ,/f + £ k 

^NO x ,/c+1 = + Q k 


M represents a transition matrix; s and Q are process errors which relate to errors in 
modeling processes, and are difficult to estimate. Since we assume the bias between 
modeled and measured N0 2 is mostly from errors in NO x emissions (Prinn, 2000; 
Napelenok et al., 2008), s and Q were set to zero. 

CAMx-DDM (Koo et al., 2007) calculates a semi-normalized N0 2 sensitivity to NO x 
emissions (unitless), as shown in Eq. (9), replacing sensitivity elements in S in Eq. (4), 


S 


N0 2 toNO x 



<?C N o 2 

^NO x 



<?C N o 2 

d((1 +*)£no x ) 


<?C N o 2 

d(1 +x) 


dC N02 

dx 


(9) 


where E represents the unperturbed NO x emission field; x is the perturbation factor. 
Hence, in this study, the DKF inversion actually seeks the optimal perturbation factor 
(x) at each iteration. The inversion processes will repeat iteratively until the perturbation 
factor for each emission region converges within a prescribed criterion, d (Fig. 2), for 
which the value of 0.01 was chosen in this study. 


2.5 NO 2 observations 


2.5.1 Satellite NO 2 measurements 


The Dutch-Finnish Ozone Monitoring Instrument (OMI) aboard NASA’s EOS Aura 
satellite, launched on 15 July 2004, is a nadir-viewing UV-Vis spectrometer that mea- 
sures solar backscattered irradiance in the range of 270 nm to 500 nm. It has been 
utilized to retrieve atmospheric N0 2 in the spectral range from 405 nm to 465 nm with 
spatial resolution down to scales of 13 km x24km at nadir view point (Levelt et al., 
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2006a, b). The EOS Aura satellite follows a Sun-synchronous polar orbit at approxi- 
mately 705 km altitude with local equator crossing time around 13:40 LT (Levelt et al., 
2006b; Boersma et al., 2007). In this study, the NASA standard product, version 2 (Buc- 
sela et al., 2013) retrieval of OMI N0 2 , gridded at 0.1° x 0.1° resolution, was obtained 
from NASA Goddard Space Flight Center and mapped to the 12 km CAMx modeling 
domain. OMI pixels with cloud radiance fraction greater than 0.5 and sizes of more than 
20 km x 63 km at swath edges were excluded in the dataset. The OMI averaging kernels 
(Eskes and Boersma, 2003) were interpolated into each CAMx model layer and then 
applied to the modeled N0 2 column density (Eq. 10), to account for the influence of 
the a priori N0 2 vertical profile used in the OMI retrieval and the OMI measurement 
sensitivities at each altitude: 


^modeled 

'NO, 


= 'Za,-x i 


( 10 ) 


where A, is the averaging kernel at pressure level /, and X ) is the CAMx modeled partial 
N0 2 subcolumn density at the corresponding pressure level. 

In order to reduce the OMI measurement uncertainties and effects from invalid data 
points, monthly averaged OMI N0 2 column densities were used in the DKF inversions. 


2.5.2 Ground and other NO 2 measurements 


The US EPA Air Quality System (AQS) N0 2 ground monitoring network data (Fig. 1) 
(http://www.epa.gov/ttn/airs/airsaqs/) were also used for inverse modeling. AQS mon- 
itors are equipped with a heated molybdenum catalytic converter that first transforms 
N0 2 to NO, and then measures the resultant NO using a chemiluminescence analyzer. 
N0 2 is then calculated by subtracting NO measured in a separate NO mode from 
the resultant NO (US EPA, 1975). Studies (US EPA, 1975; Demerjian, 2000; Lam- 
sal et al., 2008) indicate that the catalytic converter also converts fractions of other 
reactive nitrogen species (e.g. HNO s , PAN) into NO during this measurement. There- 
fore, correction factors computed from CAMx modeled concentrations by the method 
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of Lamsal et al. (2008) (Eq. 1 1) are applied before deploying the AQS N0 2 data in the 
DKF inversion: 


CF = 


N0 2 

N0 2 + JAN + (0.95PAN) + (0.35HNO 3 ) 


( 11 ) 


In Eq. (11), JAN represents the sum of all alkyl nitrates and PAN is peroxyacetyl 
5 nitrate. The CAMx model with CB05 mechanism does not output alkyl nitrates specif- 
ically, so the difference between modeled total organic nitrates and PAN was used to 
represent modeled alkyl nitrates. 

The NOAA P-3 aircraft N0 2 data (http://www.esrl.noaa.gov/csd/tropchem/ 
2006TexAQS/) and the Texas Radical and Aerosol Measurement Program (TRAMP) 
io N0 2 data, measured at Moody Tower (Fig. 1), (http://geossun2.geosc.uh.edu/web/ 
blefer/TRAMP/Final%20data/) were used to evaluate the inverse modeling results. 
The Moody Tower measurement site located at the University of Houston campus is 
approximately 70 m above the ground (Luke et al., 2010), corresponding to the CAMx 
modeling layer 2, with hourly N0 2 data available for the whole August-September 
is episode, but no coverage for the June episode. The P-3 aircraft measurement was 
made from ground level to around 5000 m height with 1 s resolution, but only available 
on 4 days (31 August, 11 September, 13 September, and 15 September 2006) during 
our modeling period. Hourly averaged aircraft N0 2 data was used to compare with the 
hourly modeled N0 2 at corresponding grid cells. Both P-3 aircraft and Moody Tower 
20 N0 2 measurements were made by using a photolytic converter, and hence did not 
require corrections via Eq. (11). 
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3 Results and discussion 

3.1 Pseudodata test for the DKF inversion with CAMx-DDM 

To evaluate the performance of the DKF inversion technique, a controlled pseudodata 
test was performed for 10 modeling days (31 May to 9 June, and 13 August to 22 Au- 
5 gust) for each modeling episode. The 10 day averaged modeled N0 2 columns at 1- 
2 p.m. from the base case were used as pseudo-observations, and the model was re- 
run with NO x emissions from each region perturbed by known factors ranging from 
0.5 to 2.0 (Table 1). Appling the DKF inversion successfully adjusted the perturbed 
NO x emissions from each region back to their base values, converging in 4 iterations 
io (Fig. 3). The robustness of the DKF inversion was tested by varying the uncertainty 
parameters, which were set to 2.0 for emissions and 0.3 for observations in the initial 
pseudodata test. While higher levels of the emission uncertainty parameter and lower 
levels of the observation uncertainty parameter led to more rapid adjustments, the final 
results of the DKF inversion were insensitive to the assumed uncertainty parameters, 
is and also to the off-diagonal elements in the error covariance matrix. 

3.2 Additional NO x emissions 

Since DKF inversions scale emissions from their original levels, an appropriate a priori 
NO x emission inventory is essential for obtaining reasonable results. The NASA Inter- 
continental Chemical Transport Experiment (INTEX-A) air quality study (Singh et al., 
20 2006) found large discrepancies between aircraft measurements and CMAQ simula- 
tions of N0 2 concentrations in the upper troposphere. Possible explanations could be 
upper tropospheric NO x sources such as lightning and aircraft NO x emissions that 
are often neglected in emission inventories. Missing NO x sources in the upper tropo- 
sphere may bias the inversion on the remaining emissions (Napelenok et al., 2008). At 
25 ground level, Hudman et al. (2010) found that the soil NO x emissions estimated by the 
widely used Yienger and Levy method (Yienger and Levy, 1995) were underestimated 
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by a factor of 2 over the United States. Therefore, in this study, the lightning and aircraft 
NO x emissions were added in the upper troposphere as described in the Sect. 2.2, and 
the soil NO x emissions were doubled from base case levels. The emission inventory 
with added lightning and aircraft NO x , and doubled soil NO x (hereafter referred to as 
the a priori emission inventory) was used for the following inversion studies. Inclusion 
of these NO x sources improves the performance of the model in simulating satellite 
observed N0 2 column densities, especially in the rural areas (Figs. 4c and 5c), and 
reduces the bias and error by around 15% (Table 3). 

3.3 Top-down NO x emissions using OMI NO 2 

3.3.1 DS inversion 

The DS inversion method was performed with OMI N0 2 column densities to create top- 
down NO x emissions for the 12 km modeling domain. The monthly averaged (3 June 
to 1 July, and 16 August to 15 September) N0 2 column densities at 1-2 p.m. were 
used to calculate the ratio of OMI to CAMx (Eq. 1). The first three modeling days 
were discarded for both modeling episodes to avoid the influence of initial conditions. 
The monthly 24 h averaged NO x emissions and modeled N0 2 column densities were 
used to determine the value of the smoothing parameter, k. In this case, k equals to 
2.0 for both episodes, indicating large influence of NO x emissions transported from 
surrounding grid cells. 

Results (Table 2) show the DS inversion scales up the NO x emissions in all seven 
regions, leading to higher estimates of modeled N0 2 column densities (Figs. 4d and 
5d) in most of the domain. However, especially in urban areas, the simulated N0 2 col- 
umn densities with inverted NO x emissions overshoot those observed by OMI. This 
indicates that the ability of NO x to influence its own lifetime via changes in OH radical 
concentrations results in significant nonlinearity between N0 2 concentration and NO x 
emission that are neglected by the DS method. Use of inverted NO x emissions does re- 
duce bias and error in simulating OMI observed column densities, while R 2 gets worse 
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(Table 3), indicating no improvement in the spatial distribution. The comparisons with 
AQS ground measurements (Table 4) indicate that the inverted NO x emissions actually 
deteriorate the simulations of ground-level N0 2 , with bias and error increasing by 70 %. 
Similar results are shown in evaluating model performance against Moody Tower and 
P-3 measurements: the DS inversion increases bias and error by approximately 30% 
and 20 %, respectively. 

3.3.2 DKF inversion 

DKF inversion using the OMI N0 2 measurements was conducted to constrain NO x 
emissions from the seven designated regions. The monthly averaged (3 June to 1 July, 
and 16 August to 15 September) OMI and CAMx N0 2 column densities at 1-2 p.m. 
were used in the inversion. All modeling grids in the inversion area were covered by 
the OMI N0 2 measurement data. The DKF inversions were performed with 2116 data 
points in one time step (1-2 p.m.). The scaling factors generated by inversion for each 
region were applied to the NO x emission inventory hourly, since we assume that the 1- 
2 p.m. N0 2 column density is contributed by the NO x emissions from all previous hours, 
and the uncertainty in the bottom-up NO x emission inventory should be the same for 
every time step. The satellite-based DKF inversions scale a priori NO x emissions by 
factors ranging from 1.02 to 1.84 in almost all regions in both episodes (Table 2), ad- 
hering to the specified uncertainty range of 0.5 to 2.0. The scaling factors tend to be 
larger over the rural and small urban regions than over the urban DFW and HGB ozone 
nonattainment regions, where the inversions scale up emissions only slightly (factors 
of 1.02 to 1.14). It results from the inversion attempts to compensate for the large 
gap between higher observed than modeled N0 2 over rural regions, despite varied 
patterns over urban grid cells. One exception occurs in the NE Texas region in the 
August-September episode (Table 2), which shows downward scaling (factor of 0.56). 
This reflects the inversion shifting emissions between NE Texas and the much larger 
surrounding N rural region (Fig. 1); taken together, the net scaling factor for the two 
regions in August-September is 1 .72, consistent with the upward scaling of rural emis- 
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sions throughout the two episodes. Apart from this anomaly, scaling factors for most 
regions were consistent across the two episodes, varying by less than 15%. 

CAMx modeled N0 2 column densities with the inverted NO x emissions (Figs. 4e and 
5e) are increased by 3-55% in all regions, but the increments are much more moder- 
ate compared to the DS method inversion. The statistical results (Table 3) indicate that 
the DKF inversed N0 2 are closer to OMI observations than the a priori case in terms 
of 20% less in bias and 10% less in error, but without improvements in the spatial 
distribution. The DS method scales up NO x emissions more than the DKF inversion 
(Table 2), making the inversed N0 2 concentrations have slightly less bias and error 
(Table 3). However, the DKF inversed N0 2 has better R 2 than that of inversed by the 
DS method, indicating the DKF inversion method has better ability to retain the spatial 
structure of NO x emissions. Each of the inversions using OMI N0 2 data actually wors- 
ens the model performance in simulating ground level N0 2 concentrations (Table 4), 
since the modeled ground N0 2 using the base case emission inventory had already 
been overestimated (Fig. 6). However, the DKF inversed N0 2 only increases the bias 
and error by 30% in simulating AQS N0 2 (Table 4), by 5% in simulating P-3 N0 2 
(Table 5), and by 8% in simulating Moody Tower N0 2 , whereas greater deterioration 
resulted from the DS inversion. 

3.4 Top-down NO x emissions using ground AQS N0 2 

Ground-level AQS N0 2 measurements were also used to drive DKF inversions of NO x 
emissions for the two modeling episodes. There are 37 ground measurement sites 
in the designated inversion regions (Fig. 1), mostly located in the urban cores. The 
N rural and S rural regions were excluded in this case because they contain too few 
measurement sites. Correction factors from Eq. (11) were applied to the ground N0 2 
before using the data in the inversion. 

The base case simulations strongly overpredicted observed N0 2 in the early morn- 
ing and late afternoon during both modeling episodes (Fig. 6), when the model may 
underestimate planetary boundary layer (PBL) heights (Kolling et al., 2013). To allevi- 

17494 


C/5 

O 

c 

C/5 

C/5 

O' 

o 

“0 

0) 

T3 

CD 


C/5 

O 

c 

C/5 

C/5 

o' 

o 

“0 

0) 

T3 

CD 


C/5 

O 

c 

C/5 

C/5 

o' 

O 

Tl 

05 

T3 

CD 


ACPD 

13, 17479-17517, 2013 

Inverse modeling of 
Texas NO x emissions 

W. Tang et al. 


Title Page 


Abstract ■ Introduction 


Conclusions | References 


Tables 


Figures 



Close 


Full Screen / Esc 


C/5 

O 

c 

C/5 

C/5 

O' 

o 

Tl 

03 

■o 

CD 


Printer-friendly Version 


Interactive Discussion 


@ © 















ate the influence from PBL heights, daily 24 h averaged N0 2 levels were used in the 
inversions. 

To address the overprediction of ground-level N0 2 , the ground-based inversions 
sharply reduce a priori NO x emissions by applying scaling factors of 0.30 to 0.57 
(Table 2). The reductions in NO x emissions reduce model error relative to the AQS 
(Table 4) and Moody Tower N0 2 observations on an hourly basis, as well as N0 2 ob- 
served by the P-3 aircraft (Table 5), but may be too sharp, as they lead negative bias 
in predicting N0 2 from the AQS monitors (Table 4) and the P-3 aircraft measurements 
(Table 5). More moderate scaling factors are obtained if the inversion is conducted with 
data only from a midday window (9 a.m.-2 p.m.) when PBL heights are less problematic 
(not shown). However, scaling factors still remain far below 1 .0 and show up to factor 
of two inconsistencies between the two episodes. 

3.5 Impacts on O 3 simulations 

0 3 concentrations and their sensitivities to changes in emissions are calculated for both 
modeling episodes using the a priori and each of the a posteriori emission inventories. 
The scaled up NO x emissions from the satellite-based DKF inversion (Table 2) lead 
to 1-7ppb higher modeled 8h (10a.m.-6 p.m.) 0 3 concentrations over most of the 
domain in the June episode (Fig. 7, top row). Largest increases occur over NE Texas 
and N rural regions (Fig. 1 ), where the a priori simulation shows 0 3 to be most sensitive 
to NO x (Fig. 7, middle row) and where the satellite-based DKF inversion scaled up 
emissions by large amounts. 

The a priori simulation shows 0 3 to be primarily sensitive to NO x over most of the 
domain, but VOC-limited in the core of the Houston region and with joint sensitivity to 
NO x and VOC in Dallas, Austin, and San Antonio (Fig. 7, left column). The satellite- 
based inversion increases NO x emissions and thus shifts the 0 3 formation chemistry 
toward being more VOC sensitive (Fig. 7, middle column). Over much of the domain, 
0 3 sensitivity to VOC increases by a factor of about 1 .5. The slight increases in 0 3 sen- 
sitivity to NO x occur because the semi-normalized sensitivity coefficients represent the 
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local slope of 0 3 -emissions response scaled to a 100% change in emissions. As the 
satellite-based inversion scales up NO x emissions, these semi-normalized coefficients 
increase, even though the impacts per ton of NO x decrease. 

The ground-based DKF inversion leads to 0 3 reductions of 3-8 ppb over urban re- 
5 gions (Fig. 7, top right), where it scales down emissions (Table 2), and less changes 
over rural regions where emissions were left unchanged due to lack of N0 2 monitors. 
The reduction in urban NO x makes 0 3 less sensitive to VOC emissions as expected 
(Fig. 7, bottom right). However, the impact on sensitivity to NO x is mixed. In urban areas 
which are transitional between NO x -limited and NO x -saturated conditions, the reduc- 
io tion in NO x emissions pushes the chemistry toward more NO x -limited conditions and 
thus increases the sensitivities. In downwind regions which are already NO x -limited, 
the sensitivities decline because there are now less NO x emissions contributing to the 
semi-normalized coefficients. 

Model performance in simulating hourly AQS ground-level observations of 0 3 indi- 
i 5 cates that the bias and error slightly worsened when each of the a posteriori inventories 
are used in place of the a priori inventory (Table 6 ). The largest deterioration comes 
from the DS inversion as the bias and error increase by around 10%, likely because 
this inversion method does not retain the spatial structure of emissions from the a pri- 
ori inventory. For the other inversions, the changes in bias and error are too slight to 
20 determine if performance is meaningfully impacted. 


4 Conclusions 

Inverse modeling has been performed using either N0 2 column densities observed by 
OMI satellite or ground-level N0 2 concentrations observed by AQS monitors to con- 
strain the NO x emissions for two regulatory attainment modeling episodes in Texas. 
25 Two inversion methods, DS and DKF, are applied to the OMI N0 2 data, and the DKF 
method is also applied to the ground-level N0 2 data. Pseudodata test results validate 
that the DKF method effectively captures known perturbations in CAMx simulations. 
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Two missing NO x sources in the upper troposphere, lightning and aircraft NO x emis- 
sions, are added into the base case NO x emission inventory, contributing 14% and 
6 % to the total NO x emissions for the June episode, and 7% and 6 % for the August- 
September episode, respectively. The underestimated soil NO x emissions are doubled 
5 from the base case, adding an additional 8 % NO x emission to the base case for both 
episodes. The additional NO x emissions increase the modeled N0 2 column densities 
mostly at rural areas and improve the inversion performance with the OMI N0 2 , but not 
with the ground N0 2 . 

The DS inversions tend to overshoot the OMI-observed N0 2 column densities since 
io this linear inversion method ignores the nonlinear influence of NO x on its own lifetime. 
The iterative approach of the DKF inversion avoids this problem, but fails to substan- 
tially improve the spatial correlation of modeled and observed N0 2 levels since it ap- 
plies only a single scaling factor to each inversion region. 

The overall tendency of the model to underpredict OMI observed N0 2 column den- 
is sities and to overpredict AQS observed ground N0 2 concentrations leads to conflicting 
results between the inversions. It is difficult to determine which observations provide 
a more reliable basis for the inversions, since none of the inversions improve model 
performance against independent data such as aircraft-observed N0 2 or ground-level 
0 3 concentrations. Whether this indicates that the a priori inventory is the best avail- 
20 able representation of NO x emissions, or that tuning of the base model led to its better 
performance, is impossible to determine. Nevertheless, this suggests that inverse mod- 
eling of NO x emissions should for now remain a complement to SIP modeling efforts 
rather than a substitute for traditional bottom-up inventories. 

The AQS ground N0 2 measurements face limitations due to the inaccuracies of the 
25 molybdenum converter method, and because the mostly urban measurement sites may 
be unrepresentative of the entire region. In addition, model shortcomings in simulating 
PBL heights may strongly bias the inversions based on ground-level observations. 

For the satellite data, several factors could explain the more spatially smeared and 
higher rural N0 2 in the satellite observations than the base model which drove the up- 
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ward scaling of emissions. Our inclusion of lightning and aircraft NO x emissions and 
doubling of soil NO x emissions narrowed but did not eliminate the discrepancy. A higher 
resolution OMI N0 2 product (retrieved with small pixels and high resolution a priori pro- 
file) has been shown to enhance N0 2 column densities in urban areas and reduce them 
5 in rural areas (Russell et al., 2011), which would more closely resemble the modeled 
distribution. Lin et al. (2012) highlighted several uncertain model parameterizations 
that impact model predictions of N0 2 column density for a given emissions inventory. 
For example, lowering the rate constant of the N0 2 + OH reaction to match the rate of 
Mollner et al. (2010) would lead to a longer NO x lifetime and reduce the gap between 
io modeled urban and rural N0 2 concentrations. Henderson et al. (2011) suggested that 
better representation of acetone and organic nitrates in the CB05 mechanism could 
help address its underprediction of N0 2 in the remote upper troposphere. Future work 
could explore how combinations of these adjustments influence satellite-based inver- 
sions. 

is The upcoming DISCOVER-AQ campaign by NASA in fall 2013 will provide verti- 
cally resolved measurements of NO x from repeated aircraft spirals in the Houston re- 
gion. This may help resolve some of the discrepancies noted here between inversions 
driven by ground-based and satellite-based N0 2 observations. The future Tropospheric 
Emissions: Monitoring of Pollution (TEMPO) mission, using a geostationary satellite 
20 with high spatial and temporal measurement capabilities, could provide a richer data 
source to drive the NO x inversions. Future work could also conduct inversions based 
on emission categories rather than emission regions, to explore potential errors in the 
emission inventory on a component rather than location basis. 
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NNX10AO05G and by the NASA Air Quality Applied Science Team. The authors thank Jim 
McKay and Doug Boyer at TCEQ for providing emission inputs, Gary Wilson and Greg Yarwood 
at ENVIRON for CAMx support, and Tom Ryerson and Winston Luke at NOAA for the P-3 
aircraft N0 2 and the Moody Tower N0 2 measurement data. 
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Table 1 . Arbitrary perturbation factors for pseudodata testing. 


Source Region 

Perturbation factor 

HGB 

1.8 

DFW 

0.6 

BPA 

1.6 

NE Texas 

0.7 

Austin and San Antonio 

1.4 

N rural 

1.5 

S rural 

0.8 
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Table 2. Scaling factors for each region from different inversions. 


Source 

Region 



3 Jun to 1 Jul 2006 




16 Aug to 15 Sep 2006 



Base NO x 
emission 
(td' 1 ) 

Priori NO x 
emission 3 
(td" 1 ) 

Scaling factor relative to priori 
(unitless) 

Base NO x 
emission 
(td- 1 ) 

Priori NO x 
emission 
(td- 1 ) 

Scaling factor relative to priori 
(unitless) 




Posteriori 
OMI-based 
DS inversion 

Posteriori 
OMI-based 
DKF inversion 

Posteriori 
Ground-based 
DKF inversion 11 



Posteriori 
OMI-based 
DS inversion 

Posteriori 
OMI-based 
DKF inversion 

Posteriori 
Ground-based 
DKF inversion 

HGB 

374 

455 

1.46 

1.12 

0.37 

382 

436 

1.56 

1.03 

0.57 

DFW 

335 

435 

1.36 

1.02 

0.34 

314 

412 

1.47 

1.14 

0.47 

BPA 

81 

97 

2.23 

1.83 

0.49 

86 

98 

2.02 

1.75 

0.42 

NE 

141 

164 

2.33 

1.84 

0.49 

155 

174 

1.69 

0.56 

0.51 

Texas 

Austin 

252 

319 

1.55 

1.28 

0.30 

248 

302 

1.82 

1.70 

0.40 

and San 
Antonio 
N rural 

522 

823 

1.93 

1.67 


543 

759 

2.00 

1.98 


S rural 

472 

728 

1.83 

1.52 

- 

489 

668 

2.04 

1.72 



a Adds lightning and aircraft NO x , and doubled soil NO x emissions to the base case. 
b Conducted with 24 h averaged ground-level N0 2 data. 
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Table 3. Performance of CAMx in simulating O Ml -observed N0 2 column densities. 


Statistical 

Parameters 


3 Jun to 1 Jul 2006 


16 Aug to 15 Sep 2006 

Base 

case 

Priori 0 

Posteriori 
OMI-based 
DS inversion 

Posteriori 
OMI-based 
DKF inversion 

Base 

case 

Priori 

Posteriori 
OMI-based 
DS inversion 

Posteriori 
OMI-based 
DKF inversion 

R 2 

0.62 

0.61 

0.42 

0.54 

0.63 

0.48 

0.40 

0.51 

NMB a 

-0.47 

-0.30 

0.087 

-0.12 

-0.54 

-0.33 

0.13 

-0.12 

NME b 

0.48 

0.32 

0.22 

0.23 

0.55 

0.39 

0.32 

0.28 


a Normalized mean bias. 
b Normalized mean error. 

c Adds lightning and aircraft NO x , and doubled soil NO x emissions to the base case. 
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Table 4. Performance of CAMx in simulating AQS Ground-level N0 2 a . 


Statistical 3 Jun to 1 Jul 2006 1 6 Aug to 1 5 Sep 2006 

Parameters 



Base 

case 

Priori 

Posteriori 
OMI-based 
DS inversion 

Posteriori 
OMI-based 
DKF inversion 

Posteriori 
Ground-based 
DKF inversion 

Base 

case 

Priori 

Posteriori 
OMI-based 
DS inversion 

Posteriori 
OMI-based 
DKF inversion 

Posteriori 
Ground-based 
DKF inversion 

R 2 

0.55 

0.54 

0.49 

0.51 

0.52 

0.50 

0.49 

0.43 

0.44 

0.47 

NMB 

0.82 

0.91 

1.66 

1.30 

-0.16 

0.35 

0.42 

1.10 

0.72 

-0.24 

NME 

0.96 

1.03 

1.71 

1.37 

0.48 

0.63 

0.67 

1.21 

0.91 

0.48 


a Hourly AQS data was used to compare with modeled N0 2 at corresponding locations. 
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Table 5. Performance of CAMx in simulating P-3 aircraft-observed N0 2 . 


Statistical 1 6 Aug to 1 5 Sep 2006 a 

Parameters 



Base 

case 

Priori 

Posteriori 
OMI-based 
DS inversion 

Posteriori 
OMI-based 
DKF inversion 

Posteriori 
Ground-based 
DKF inversion 

R 2 

0.23 

0.23 

0.25 

0.22 

0.22 

NMB 

0.10 

0.10 

0.46 

0.15 

-0.14 

NME 

0.99 

0.99 

1.24 

1.02 

0.86 


a Comparison available for only four days (31 August, 1 1 September, 13 September, and 
15 September 2006). 
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Table 6. Performance of CAMx in simulating AQS hourly ground-level 0 3 . 


Statistical 3 Jun to 1 Jul 2006 16 Aug to 15 Sep 2006 

Parameters 



Priori 

Posteriori 
OMI-based 
DS inversion 

Posteriori 
OMI-based 
DKF inversion 

Posteriori 
Ground-based 
DKF inversion 

Priori 

Posteriori 
OMI-based 
DS inversion 

Posteriori 
OMI-based 
DKF inversion 

Posteriori 
Ground-based 
DKF inversion 

R 2 

0.61 

0.61 

0.63 

0.58 

0.50 

0.52 

0.51 

0.46 

NMB 

0.01 

0.12 

0.02 

0.04 

0.38 

0.49 

0.41 

0.40 

NME 

0.29 

0.37 

0.30 

0.30 

0.47 

0.58 

0.50 

0.48 
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Fig. 1. 12 km CAMx modeling domain for eastern Texas (black square), inversion regions 
(shaded), ground AQS N0 2 monitoring sites (blue triangles), and Moody Tower (red circle). 


g 

CO 

o 

c 

co 

co 

O' 

o 

"0 

O 

T3 

CD 


D 

CO' 

o 

a 

CO 

c o 
o' 

13 

“0 

o 

T3 

CD 


D 

cn' 

o 

c 

CO 

CO 

o' 

o 

"0 

o 

T3 

CD 


g 

CO 

o 

c 

CO 

CO 

o' 

o 

"0 

a> 

■o 

CD 


Inverse modeling of 
Texas NO x emissions 

W. Tang et al. 


Title Page 


Abstract 

Conclusions 

Tables 



Back 


Introduction 

References 

Figures 



Close 


Full Screen / Esc 


Printer-friendly Version 
Interactive Discussion 



17511 















ACPD 

13, 17479-17517, 2013 



Fig. 2. Schematic diagram of Kalman filter inversion process. 
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a) b) 



Base modeled N0 2 (10 5 molecules/cm 2 ) Base modeled N0 2 (I0 15 molecules/cm 2 ) 

Fig. 3. Pseudodata test showing that the DKF inversion accurately adjusts the NO x emissions 
from the perturbed case (a) to the a posteriori case (b) to match the desired base N0 2 column 
densities. Similar performance is found for the 13-22 August test period. 
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a) 


b) 


c) 



d) e) 



Min- 0.747 at(l, 18), Max- 17.517 at (37,34) Min- 0.611 at (1.14). Max- 10.223 at(33,10) 


Fig. 4. Monthly averaged (3 June to 1 July) tropospheric N0 2 vertical columns at 1-2 p.m. from 
(a) OMI observations, and from CAMx simulations using (b) base case emissions inventory, 
(c) a priori emission inventory (with additional lightning, aircraft, and soil NO x ), and OMI-based 
inverted NO x emissions using (d) DS and (e) DKF methods. 
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a) 


b) 



Min= 0.913 at (37.1), Max= 6.308 at (18,38) 


46 


1 

Min= 


46 

0.375 at (1,15), Max= 8.617 at (33.10) 



d) e) 



Min= 


0.507 at (1,15), Max= 9.461 at (33,10) 


46 

13.520 at (18,17) 


6.000 46 

5.143 

4.286 

3.429 

2.571 

1.714 

0.857 


^ 0.000 ! 
lO^ISmolecul^s/cm^ 


at (1,17), Max= 



Fig. 5. Same as Fig. 4, but for the August-September episode 
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c) 



Mill- 0.434 at(1.15J, Max- 8.944 at(33,10) 
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Time (hr) 

Fig. 6. Daily variations of modeled (solid line) and observed (dashed line) ground N0 2 concen- 
trations for the June (red) and August-September (blue) episodes. Note: N0 2 concentrations 
were taken by averaging monthly data for all sites. 
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a priori 


u postmori:OMI bused UK* inversion 


a posteriori: Ground based DKf inversion 



0.073 16 

0.009 

0.066 

0.062 

0.059 

0JJ55 

0.05? 

0.000 1 
ppmV 


1 <16 
Min- 0.0<W at (1 .1 \ Max- 0.0/3 at (13,H) 


i- -0.000 at (1 6.40). Max- 


-16 

-0.001 at (46.32) 



Min- -0.000 at (46.42). Max- 0.012 at (33.10) Min- 0.000 at (<16.73). Max- o.ooi at (33.36) Min- -0.00/ at (33.10). Max 0.000 at (16.31 ) 


Fig. 7. Monthly 8h (10a.m.-6p.m.) averaged ground-level 0 3 concentrations (top), 0 3 sensi- 
tivity to NO x (middle), and 0 3 sensitivity to VOC (bottom) for the a priori case (left column), and 
differences (a posteriori minus a priori) for the OMI-based (middle column) and ground-based 
(right column) DKF inversions in the June episode. The August-September episode shows 
similar results. 
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